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Abstract 

Demonstrating the completeness of wave functions solutions of the radial Schrodinger equation is a very 
difficult task. Existing proofs, relying on operator theory, are often very abstract and far from intuitive com- 
prehension. However, it is possible to obtain rigorous proofs amenable to physical insight, if one restricts the 
considered class of Schrodinger potentials. One can mention in particular unbounded potentials yielding a 
purely discrete spectrum and short-range potentials. However, those possessing a Coulomb tail, very important 
for physical applications, have remained problematic due to their long-range character. The method proposed 
in this paper allows to treat them correctly, provided the non-Coulomb part of potentials vanishes after a 
finite radius. Non-locality of potentials can also be handled. The main idea in the proposed demonstration 
is that regular solutions behave like sine/cosine functions for large momenta, so that their expansions verify 
Fourier transform properties. The highly singular point at k = of long-range potentials is dealt with properly 
using analytical properties of Coulomb wave functions. Lebesgue measure theory is avoided, rendering the 
demonstration clear from a physical point of view. 

1 Introduction 

Completeness of wave functions of the Schrodinger equation has always been a difficult subject in quantum 
mechanics. Due to the importance of eigenfunctions expansions in quantum physics, a global comprehension of 
the phenomena underlying the property of completeness is necessary. 

Demonstrations of completeness have been treated by many authors. Fourier and Fourier-Bessel series 
come back to the works of Dirichlet for the former [TJ and Hankel, Schlafii and Young for the latter [2] (see 
Ref. [3] for a thorough study of Bessel functions properties and demonstration of completeness). Both proofs 
use complex analysis only. Fourier series demonstration is standard and rely on Dirichlet kernel properties. 
The Fourier-Bessel series proof is similar to the former through the use of a generalization of the Dirichlet 
kernel. The general case of differential equations defined in a finite interval, of which Fourier and Fourier- 
Bessel expansions are a particular case, is treated within regular Sturm-Liouville theory [H[S]. Completeness 
can therein be demonstrated using Rayleigh's quotient [6j [7] or closeness arguments of the considered set of 
eigenstates with another basis [4], with which, however, point convergence cannot be handled. The regular 
Sturm-Liouville problem with non-local potentials has been treated by Fubini and Lichtenstein [8]. Extension 
to the singular case of infinite intervals, albeit for differential equations only, has first been considered by Weyl 
[§] , and continued afterward by Titchmarsh [5] . The study of spectral theorem of the general linear operator can 
be found in Ref. [10] . If one restricts operators to be short-range spherical potentials for the radial Schrodinger 
equation, the method of Newton [11] can be applied. It deals with complex integration only and, thus, has 
the advantage to be straightforward mathematically. Its extension to the Coulomb case was, however, not 
considered. It has been effected recently, but only for the case of Coulomb wave functions, for which analytical 
properties of the confluent hypergeometric function can be employed (T3] • 

Newton's completeness relation has played a fundamental role in the development of non-hermitian formal- 
ism, as it is the starting point of the Berggren completeness relation [14] , where bound, resonant and complex 
scattering states are used. It has been shown to be a powerful tool to expand strongly correlated nuclear states 
[15 ] [TB I IT7] . With a continuous spectrum, the theory of rigged Hilbert spaces [18, 19] is necessary to be able to 
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use Dirac formalism correctly. This theory cannot be built unless completeness of eigenstates wave functions 
has been proved. 

Newton's completeness relation extended to the proton case with repulsive Coulomb potential was demon- 
strated in Ref. [T7] , using discrete completeness relations given by box boundary conditions, becoming continuous 
letting the radius of the box go to infinity. However, mathematical details were missing, which could not be 
treated in Ref. [17] as the principal motive of the paper was the study of loosely bound and resonant nuclear 
systems and not the completeness demonstration by itself. They will be described in this paper. The aim of 
this paper is also to extend the demonstration of Ref.[T7] to potentials with attractive Coulomb tails. The 
main difficulty therein is to treat properly the infinite set of bound states accumulating at zero energy for 
attractive potentials. The generalization to the non-hermitian case of complex potentials will be also be per- 
formed. Standard mathematical results of fundamental importance for our demonstration will be presented in 
appendices. 

2 Hamiltonian potentials 

One will analyze the sets of eigenstates provided by one-body Schrodinger equation Hamiltonians. One will 
call eigenstates both bound and scattering states for simplicity, even though the denomination is improper for 
non-square integrable positive energy states. We will first consider the radial Schrodinger equation for real 
potentials, non-hermitian complex potentials being handled afterward. Potentials can be non-local, with their 
local part denoted by v(r) and their non-local part by w(r,r'). They are demanded to vanish identically for 
r > Rq, r' > Rq, except for pure Coulombic asymptotic, which can be attractive or repulsive. v(r) and w(r, r') 
uphold the following conditions: 

(1) 

and w(r, r') -> , r -> R Vr' e [0 : R ], (2) 

where V c S R. v(r) will be supposed to be continuous except maybe for a finite set of radii, where it has to 
be integrable, and w(r,r') integrable on its domain of definition, unless explicitly stated. w(r,r') is naturally 
supposed to be symmetric in r and r' . The first condition in Eq.Q is necessary for eigenstates to be well-defined 
and is automatically fulfilled for Hartree-Fock potentials [5D] . Demanded conditions are virtually always verified 
in practice, as potentials subtracted from their Coulomb component usually go very quickly to zero at large 
distance. The studied Schrodinger equation then reads: 

u"(r) — ( — — -rj- — - + v(r) — k 2 J u(r) + [ w(r, r') u(r') dr' (3) 
V r J Jo 

where u(r) is the radial wave function, not necessarily eigenstate of the Hamiltonian, defined by v(r) and w(r, r'), 
I is its angular momentum and k is its linear momentum. I can bear non-integral values in our demonstration 
but is restricted for the moment to t > 0. Negative values of £, for which one can always consider £ > — i due 
to invariance of Eq.© with I — > — I — 1, will be effected in the context of non-hermitian potentials. Existence 
and unicity of u(r) functions in Eq.([3]) for arbitrary k is provided for r > Rq by the Cauchy-Lipshitz theorem, 
as there Eq.Q is local and v(r) continuous except at mentioned singularities, whose treatment is standard. For 
r < Rq, it can be equivalently written as an integral equation which then can be handled by Fredholm theory 
as Rq is finite and v(r), w(r, r') are integrable. Hamiltonians bearing an effective mass are implicitly accounted 
for with Eq. (|3|) as they can always be rewritten to verify a Schrodinger equation bearing no effective mass via a 
point-canonical transformation [21j . The considered effective mass must evidently become constant for r > Rq. 

The Hamiltonian defined in Eq.j3|) in open radial interval possesses a lower bound for eigenstate energy, 
i.e. a bound ground state or a continuum of positive energy states only. This property is immediate for local 
potentials, with which convexity of wave functions for sufficiently large —k 2 values is used to demonstrate this 
property [7j- As non-locality in our class of potentials vanishes after i?o, it is possible to use this type of 
argument for non-local potentials. Indeed, if a wave function, regular in r = but not necessarily for r — » +oo, 
has its last node in [Rq : +oo[, making its energy more negative will either send this node to infinity and have it 
disappeared or make it converge to a finite radius Ri > i?o- Hence, there exists in all possible cases a negative 
energy e m i n and radius R\ > for which all wave functions sustaining e < e m i n are nodeless in [Ri : +oo[. Thus, 



v (r) = — , r > Rq, 
r 



w(r, r') ~ w(0, r')r 



, r 
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the Hamiltonian box spectrum for R > R\ has a ground state which will converge either to the Hamiltonian 
ground state in open radial interval or become part of the continuum if the open radial interval spectrum has 
no discrete part therein. 



3 Large momentum expansion of wave functions 

The fundamental idea of the proposed demonstration is that wave functions possess an analytical expansion in 
the limit of large momentum. Formulas are based on the semi-classical approximation inspired from Wentzcl- 
Brillouin-Kramer (WKB) theory [22] • In Ref. [17] . however, spherical Bessel functions are used instead of the 
usual sine/cosine functions. Their use, indeed, allows to take into account properly the first turning point 
induced by the centrifugal barrier in Eq.Q, so that it is available for all r > 0, k sufficiently large. One will 
demonstrate its validity in this section. The simplest way for this purpose is to introduce the error equation 
[23) , i.e. one writes the exact solution of Eq.Q as a product of approximated solution and error term. We 
will suppose for the moment that v(r) and w(r, r') are twice differentiable. In order to study our semi-classical 
approximation of Eq.Q, we define the following functions: 

A(r) = yjk 2 - v{r) , A(r) = / A(r') dr', (4) 

Jo 

u fo {r) = A(r)"* //(A(r)) , u f (r) = u fo (r) [1 + e(r)}, (5) 

where k is the linear momentum of the wave function u(r) of Eq.Q, fi is a spherical Riccati-Bessel function, 
denoted as jg, hf or hj respectively regular, outgoing and incoming wave solutions of Eq.Q for v(r) and 
w(r,r') identically equal to zero, Uf (r) is the approximated solution and e(r) the error function of Eq.®. k 
can be arbitrary large (which one will denote also as k > k m i n ) so as to have A(r), A(r) and u/ (r) defined 
Vr > 0. For the method to work, one needs uj (r) never to vanish except maybe at r = 0. Hence, one will 
have fi = ji if r € [0 : r ], r chosen small enough so that /^(A(r)) ^ for r g]0 : r ], and fi^hf\ir> r . 
Inserting Eq.([5]) in Eq.([3]), one obtains the error equation verified by e(r) with r > 0: 

d 



e"(r) + e'(r) 



— \og(u h (xY) 



= F(k,r,e), (6) 



F(k, r, e) = [1 + e(r)] I -A" (r)A(r)- 1 - -A'(r) J A(r) 
+ i(i +m + (i r)](\- k2 - v{r) 



r 2 A(r^ 
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w{r,r')u fo (r')[l + e(r')} dr'. (7) 



Ufa ( r ) -'O 

Note that the term depending on £(£ + 1) appears in Eq.([7j) because of the use of Bessel functions in Eq.©. 
One can solve Eq.© as a first-order equation according to e'(r) treating F(k,r,e) as formally independent of 
e(r). This provides the integral equation verified by e(r): 



r dr' r 

< r )= -^Tlv u f0 (r") 2 F(k,r",e)dr", (8) 

Jr d u fo V ) Jr d 



where r^ — for fi — ji and r c i = +oo for fi = hf . The iterative Picard's method to formally solve Eq.(| 
comes forward: 



dr' f r 

e (r) = , e»+i(r) = / f ^ 2 / u fo (r") 2 F(k,r" ,e n ) dr" , n > 0, (9) 



d u fo(r') 



with which e' n+ i{r) is obtained by a simple differentiation. One will show that e n (r) — > e(r) for n — > +oo and 
e(r) = 0(fc -2 ), e'(r) = 0(fc -2 ) when fc — > +oo by recurrence on n G N, and that Vr > 0. £o(r) — 0(k~ 2 ) and 
e' (r) = 0(k~ 2 ) are trivially verified. One thus supposes e„(r) = 0(k~ 2 ) and e^(r) = 0(k~ 2 ) when fc — » +oo, 
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defined Vr > 0, and this property has to be demonstrated for e„ + i(r) and e' n+1 (r). F(k, r, e) = 0(log(r)r -3 ) for 
r — > +00 in Eq.([7]), so that e n +i(r) is defined with r^ = +00 in Eq.®. The non-local term involving w(r,r') 
in Eq.© converges in r' = and has a finite limit for r — > because of Eq.©, U/„(r) ~ r^ +1 or U/ (r) ~ r - ^ 
equivalents (up to one constant) for r — > and the fact that e n (0) is finite. e n +i(r) in Eq.© remains finite 
for r — > 0, even though F(Jc, r, e„) is undefined in r = by way of the term proportional to £(£ + 1), which is 
0(r _1 ) for r — ► 0. This is demonstrated using Eq.© and u/ Q (r) ~ r £+1 or Uf (r) ~ r~^ equivalents for r — > 
as well. The same argument holds for e^ +1 (r). From Eqs. (|4l7p . one can see that F(k,r,e n ) is the sum of two 
terms, one being the product of 1 4- e n (r) multiplied by a function easily shown to be 0(k~ 2 ) when k — > +00, 
and the other involving only the non-local potential w(r, r'). Due to finiteness of e n (r) for r > 0, the first term, 
integrated in Eq.©, will remain 0(k~ 2 ) when k —* +00. The other term involving w(r,r') is treated as in 
Ref. |17j. with two integrations by parts, using Eq.©: 
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w(r,r') Uf a (r') dr' = — - 



Uf (r) Jo ' «/ (r) 



dw 1 r ' =Ro 1 f R ° d 2 



where Uf {r) is defined so that Wl(r) = w/ (r), {Mf (ro),Wf (ro)} chosen so as to have Uf (r)uf (r) _1 = 0(k~ 2 ) 
for k — > +00, choice rendered possible via the use of fg asymptotic behavior for large arguments [24] . As w(r, r') 
is bounded as twice differentiable on [0 : Ro] 2 , the non-local term of Eq. (fT0f is 0(k~ 2 ) when k — * +00. The 
remaining part of the non-local term in Eq.© involving e n (r') reads after one integration by part: 



« A( r) ^ dr '=-^r^ {ry)w dr> > (U) 

where w e<l (r,r') — w(r,r') e n (r'). As U'j Q (r)u ^(r)^ 1 = 0(fc _1 ) and as it is assumed that e n (f] = 0(k~ 2 ) and 
e' n (r) = 0(k~ 2 ) when k — > +00, the integral in Eq. fTTj) is 0(k~ 3 ) when fc — > +00. Thus, e n+ i(r) = 0(k~ 2 ) and 
£n+i( r ) = 0(k~ 2 ) when A; — > +oo. Using Eqs. (l6l7l9|10|ll|) . the following equation holds using recurrence on 
n e N: 

e «W = : %^ + °(^) ,neN,fc->+oo, (12) 

where £(r, fc) is independent of n and bounded \/k > K , K large enough, by a constant independent of fc, while 
0(fc -3 ) is bounded Vrt £ N by a constant independent of n. Thus, for k > K, K large enough, e(r) can always 
be written by way of a converging Neumann-Liouville power series in k^ 1 . Hence, as all terms in Eq. (|12p can 
be majored by constants independent of n, Eq. (fT2| is valid for e(r) as well. 

Noting the clear linear independence of both u / (r) functions involving fg — for k large enough, one can 
now express u(r), eigenstate of Eq.©, in terms of the Uf(r) functions of Eq.©. u(r) has to vanish in r = 0, so 
that one can write without loss of generality: 

u(r) = M A(r)~5 \je(A(r)) + e(r)} , r < r , (13) 
«(r) = A(r)~* [Af + (hj(A(r)) + e+(r)) + AT (ftj(A(r)) + e"(r))] , r > r , (14) 

where A/", A/" + , A/" - are normalization constants and e(r), e + (r) and e _ (r) the functions defined in Eq.© for 
respectively fg = jg in Eq.(fT3"]) and fg — hf in Eq. (fT4]) . multiplied by fg for convenience. By continuity of u(r), 
its expressions in Eq. (fT3|) and Eq. (fT4f must be equal at ro > 0, which can be chosen to be arbitrarily small. 
It readily implies that N = 2iJ\f+ = -2iAf~ and 2ie(r) = e+(r) - e~(r) in Eqs. (|13I14|) . Thus, e(r) = 0(fc" 2 ) 
Vr > when k — > +cxd in Eq. (fT3|) . Eq. ([T3| provides the asymptotic expansion in fc _1 of Ref. [17]: 

2k 



u(r)=C k j i (kr)-C k -^j' e (kr) + 0(-£) , fc — > +00, (15) 



V(r) = / u(r') dr', (16) 



where Cfe is a normalization factor. u(r) and the two first terms of the right-hand side of Eq. (fT5]) are bounded 
Vr > and k — > +00 by a constant depending only on V(r) and Ck- This property immediately holds for 
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O (%t) in Eq. (fT5|) . Thus, Eq. p3|) is still valid when lifting the condition of twice differentiability for v(r) and 
w(r,r') for (r, r') € [0 : Rq] 2 , as the space of twice-differentiable functions is dense with respect to 1-norm in 
the set of potentials verifying the conditions stated in Sec.©- If fc — > +00 and r — » +00, it is preferrable to use 
the following expansion, deduced from Eq. (fT2|) and Eq. (|T3jl as well: 

u(r) = C k UHr))+C k ^P-, (17) 
where sup k>k . n \a k (r)\ — > when r — » +00. 



4 Cancellation of Dirac delta distributions 

As described in Ref. [T7] , completeness relations involving continuous set of states are handled with box discrete 
completeness relations, whose radius R will go to +00. Used box boundary conditions are u(0) = u(R) = 0, 
which define uniquely eigenstates wave functions. From the fact that the Hamiltonian of Eq.Q possesses therein 
a ground state (see Sec.©), the completeness of a discrete set of eigenstates of Eq.([3]) defined in a finite interval 
[0 : R] can be procured for scalar product norm [6l [7] . Pointwise convergence for our class of Hamiltonians will 
be a consequence of the results obtained in this section. 

The completeness relation of eigenstates of the Hamiltonian of Eq.Q with box boundary conditions in a 
finite interval [0 : R] formally reads: 

+00 

y]u n (r)u n (r f ) + ^ u ( k rn,r)u(k m ,r')(k m - k m -i) = S(r - r 1 ) (18) 

neb m=l 

where the u n {r) states are called "bound" , as they will become square-integrable on the real axis for R — » +00, 
while the u(k m ,r) states are called "scattering", as they will converge to non-integrable scattering states for 
R — > +00. Bound states of positive energy, called bound states embedded in the continuum (BSEC) [25] . 
are allowed, but bound states of energy zero are forbidden. This last point is clearly of no consequence for 
practical calculations. The number of BSEC's for our class of Hamiltonians will be demonstrated to be finite 
and independent of R as long as R > Ro in Sec.©, so that they will pose no problem when R — > +00. Wave 
functions are normalized by way of the following equalities: 

f R 1 

u n {r) 2 dr = 1 , / u(k m ,r) 2 dr = - — — , (19) 

where k m — fc m -i is used instead of k m+ \ — k m in Rcf.[17J as it is more convenient for attractive Coulomb case. 
Intuitively, for R — > +00, the sum of Eq. (fT8|) involving bound states will remain finite (repulsive/no Coulomb 
case) or become infinite (attractive Coulomb case), while its infinite series built from u(k m , r) states will become 
an integral. One can also foresee the appearance of Dirac delta normalization of scattering states in Eq. (|19p . 
The limiting process cannot be done directly, however, as Eq. (TT8"]) is purely formal due to the divergent character 
of its series. In order to avoid convergence problems, a new method has been introduced in Ref. [17] . It consists 
in subtracting the completeness relation generated by Fourier-Bessel series to the one of Eq. (fT8|) . One then 
generates a convergent series for fixed (r, r'), so that the transformation from series to integral when R — > +00 
can be effected properly. 

One defines Bessel discretized states and formal completeness relation analogously to Eqs. (|18|19|) : 

+00 

^ B 2 Km jt{K m r)ji{n m r')(K m - K m -i) = 5{r - r') (20) 

m— 1 

B L / j/M 2 dr = - — - — ■ (21) 

where B Km is a normalization constant and K m , m G N* , is a discretized linear momentum of the Fourier-Bessel 
series, verifying ji(n m R) = 0. One recalls the asymptotic expansion of jt{x) when x — > +00: 
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where at = -£(£ + 1). Subtracting Eg. ([20]) from Eq.([T8]) provides: 
S R {r,r') = ^u n (r)u n (r') 



+ 00 



+ ^2 \u(k m ,r)u(k m ,r')(k m - fc m -i) - Bl m je(K m r)jt(n m r')(K m - n m -i) 



(23) 



where S R (r, r') is for the moment a distribution only known to be equal to zero in a weak sense. It is necessary 
to show that the series of Eq. (|2"5)) converges V(r, r') G [0 : R] 2 , after which one can infer that Sn(r, r') = in a 
strong sense. For this, one will rewrite the main results of Ref. [17j. 

We first determine the asymptotic expansion of k m and K m when m — > +00. From Eqs. (|15l22[) . one obtains 
for m — * +00: 



(m + §) tt at + RV{R) 



k m = — 'J' - + - 1"' + O [ ) , K r , 



R 



IRm-n 



(m + |) 7r at 



R 



2Rmn 



O I — 

1 * 



1 

m 1 



(24) 



The condition that the same integer m enters both linear momenta expansions of Eq. (|24|) without finite shift is 
enforced demanding that u(k m , r) and jt(Kmf) have both m nodes, having for consequence that a finite number 
of u(k m , r) in Eq. (|23[) vanish if bound states are present therein. The Ck m and B Km constants expansion defined 
with Eqs. (|15ll9l2ip are calculated with Eqs. (|15I22|) . using the fact that the integral of Eq. (|2"Tj) is analytical: 



-B 2 (K m R)jt-l(K m R) 



R(n m K m — 1) 



c. 



jt(k m R) 2 - jt+i(k m R)jt-i{k m R) 
2 



R 2V(r) , 



"TO 



km. R 



jt(k m r)j' e (k m r) dr + O ( — 



R(k m fe m _x) 



(25) 



(26) 



Eq. (f2"2")) and the fact that the integral involving V(r) in Eq. (f2"6")l is 0(fc m 2 ) when k m — > +00, which can be seen 
with a partial integration, provide for m — > +00: 



(27) 



Ck m = \J^ + 0(m- 2 ) , S Km = ^ + 0(m" 

The general term of the series of Eq. ([23| reads for m — ► +00, by way of Eqs. (|15|22|24|27| : 

u{k mi r)u{k m ,r'){k m - fc m _i) - B 2 m jt(K m r)jt(K m r')(K m - K m -i) 

= ^)-v(^))-(^^) sin ^ (r _ /) + „ 



fl(V(r)+V(rO)-(r + rQV(i?) . /vrm 

sm — — (r + r J + £tt 



2mnR 



R 



2R 
r + r' 
2R 



0(m~ 



(28) 



The non-absolutely convergent part of Eq. ([28|) consists in a standard Fourier series term. Abel transformation 
can be performed on Eq. ((28|) in order to have it normally convergent with respect to r and r' (see App.([A])). 
As a consequence, Su(r, r') is finite and continuous V(r, r') G [0 : R] 2 . Thus, owing to Eqs. (|23l28[) : 

rR rR +°° I" rR 

\ S R (r,r') tt(r') dr' = y\ w (r) / u n (r')u(r') dr' + ^ u(&m, r)(k m - k m -i) / u(k m ,r')u(r') dr' 
J ° neb Jo m=i L 70 



B l m je(n m r)(K m - K m _i) / jt{n m r')u{r') dr' 



u{r) — u(r) = 0, 



(29) 



where u(r) is a state equal to u n (r), n G N or u(k m ,r), m G N* (see Eq.([TH])). This result is obtained by way 
of orthonormality of the set of it ra (?") and u{k m , r) states, pointwise completeness of the set of Bessel functions 
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and by inverting series and integral integrating Eq. (j2"3]) between and R, which is justified due to the normal 
convergence obtained by Abel transformation of the series of Eq. ([23|) . As the set of u n (r) and u(k m , r) states is 
complete for scalar product norm, J Q Sr(t, r') 2 dr' = 0. As Sr(t, r') is finite and continuous V(r, r') £ [0 : R] 2 , 
S R (r,r') = V(r,r') £ [0 : R} 2 . 

One will show that it implies pointwise completeness of the set of u n (r) and u(k m , r) states in Eq. (|18p . One 
considers an arbitrary wave function /(r), assumed to yield a Fourier-Bessel series on [0 : R\. Its expansion with 
the set of u n (r) and u{k mi r) states is equal W £ [0 : R] to its Fourier-Bessel series expansion, which can be seen 
calculating J* Q fl Sii(r, r') f(r') dr' similarly as in Eq. ([2l)]) . The set of u n (r) and u(k m , r) functions then possesses 
the same properties as Fourier-Bessel, and hence Fourier series [3J, i.e. the generalized Fourier expansion of f(r) 
provided by Eq.® is equal to lim^o f ir+s) + f{r ~ s) . 



5 Wave functions behavior for box radius R —> +00 

In order to be able to let R go to +00 in Eq. (|2"3"]) . one has to show first that wave functions converge to their 
open radial interval counterparts. The case of bound states of negative energy can be treated considering the 
appearance of nodes at r = R when \k 2 \ decreases in Eq.© [7], as Eq.© is local at large distance (see Sec.©). 
The fact that they are square- integrable on [0 : +00 [ readily implies uniform convergence of a box bound wave 
function u n (r), n £ N fixed, with respect to r £ [0 : +00 [ when R — » +00 to its limit defined in open radial 
interval. The number of bound states of negative energy will become infinite only for potentials with attractive 
Coulomb tail in our class of Hamiltonians. Indeed, it occurs only if the number of nodes of the regular solution 
u(k = 0, r) is itself infinite [22]. This can be proved with WKB approximations of u(k = 0, r) for r — ► +00 (see 
Sec. (|6.ip for the Coulomb attractive case). 

The case of scattering states is different, as one has to show both the accumulation of box scattering states 
in ]0 : +00 [ for R —> +00 and the consistency of their normalization by way of Eq. (fTT))) . For their study, one 
considers a linear momentum k > and R always chosen larger than Rq. Eq.© implies that the wave function 
u(k, r) of linear momentum k > behaves for r > Rq as: 



„( fc) r ) = ^J-Mk (s+h+(kr) + S-hJ v (kr)) , (30) 

where 77 is the Sommerfeld parameter, here equal to h lr] is the outgoing/incoming irregular Coulomb wave 
function [53] and A4, S + and S~ are normalization constants chosen so that A4 £ M and S + S~ — \ (reality of 

2 



wave functions imply that S + = (S )*). The J f factor is introduced in order to have A4 — > 1 for scattering 
states when R — » +00. Due to Eq.(|5D]). wave functions for k > are not square integrable in [0 : +oo[ unless 
■A/fe = therein. This implies that all BSEC's of the open radial interval problem are present in the box spectrum 
and independent of R once R > Rq, as they have to uphold box boundary conditions for R = Rq and vanish 
identically afterward. Moreover, it was demonstrated in Sec© that all wave functions of linear momentum 
k large enough verify Eq. (|15|) for fixed box radius i?, so that it is impossible for wave functions of arbitrarily 
large k to vanish identically in [i?o : R] and not in [0 : Rq]. Thus, BSEC's have to be in finite number and can 
be all put in the bound states sum of Eq. fTS]) . We can thus always suppose that k is not a BSEC momentum, 
implying Afk 0. Using asymptotic expansion of Coulomb wave functions for large arguments in Eq. ([3Ul) 24 , 
one obtains: 



u(k,r) = \I-Af k 



£(£ + !) ( 1 

sin(£;r - ?7log(2fcr) + 8(k)) + — cos(fcr - ?7log(2fcr) + 6(h)) + O 

2kr 



(31) 



with S(k) £ [0 : 2ir[ the phase shift associated with u(k, r) and the k dependence is written explicitly in the rest 
term to allow k — > +00 along with r — * +00 at the end of the calculation. For u(k n ,r) eigenstate of Eq.(j3j, 
n £ N*, implying u(k n ,R) = 0, the following equality holds from Eq. ([5T]) : 



_ nn-6 n V c log(2k n R) 1(1+1) ^ f 1 \ 
ICn ~ R + 2k n R 2k n R +U \k*&)> {6} 

where r\ n and S n are its respective Sommerfeld parameter and phase shift, nevertheless modified so that one 
might have 6 n ^ 8(k n ). Indeed, in Eq. (|3"2"]) . S n has been defined so that Eqs. (|24l32p provide the same value of 
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linear momentum for k n and k m if one has n = m in both equations. For R — > +00, k n — ► fc and 5 ra converge 
to a finite value for n £ N appropriately chosen, as u(k n ,r) differs significantly from its asymptotic value on 
Ea.(|3ip only in a finite region of the real axis. One can see from Eq.(|32p that for R — * +00: 



o(R), 



(33) 



as supposing therein the opposite for sup(fcj,)i{<>i? leads to a contradiction. For R large enough, it is thus 
possible to choose n so that k n < k < k n+ i. The difference fc n +i — k n reads by way of Eq. (|32p : 



Hence, k n+ i — k n 
k n — » k when R — > 



> in Eq 

-co, Eq, 



7T 



>n+l 



Sn 



R 



K_ (\og{2k n+l R) 
2R \ 



2R 



1 



O 



klR 2 



for R — ► +00, and consequently <5 n +i 
can be simplified to: 

7T / 1 



log(2fc„i?) 



5 n — > as well. As 



«„+i 



(34) 
k and 

(35) 



This proves that the set of {fc n } n6 N* becomes uniformly dense in ]0 : +oo[ when R — » +00 (see Eqs. (|33l34l35[l ). 

We can now show that Eq. (flll|) provides the standard Dirac delta normalization of scattering functions for 
R — » +00. Inserting Eqs. (|31l35[) in Eq. (JTSJ) , one can derive for R — > +00: 



Ho 



sin(fc„r - 77 log(2fc„r) + 5(fc„)) + O [ — !— 



dr + 0(l) 



= § + o(R), 



(36) 



where the part of the integral from to Rq is accounted for by the O(l) term of the left-hand side of Eg . (f3"6"|) . 
The latter integral is bounded when R — > +00, and hence k n — > k, because u(k n , 0) = and M k u{k n , Rq) has 
a limit therein from Eq. (|30p . implying A/jT 1 u(fc,r) boundedness in [0 : i?o] for -R ~~ * +co. The left-hand side 
integral in Eq. ([36"P is dealt with similarly to Eq. (f26]) . Consequently, we arrive to the standard value of A4 of 
scattering states for R — > +00, equivalent to Dirac delta normalization: 



(37) 



A4 = l + o(l), 

uniformly in all intervals [kd : +00 [, kd > 0, and in ]0 : +00 [ if I — V c = as only the principal term of Eq. (l3ip 
remains in this case. As all solutions of Eq. (|3j regular at r — for fixed k and R > Rq are proportional, one has 
proved the convergence of box scattering states to the continuum of positive energy states of Eq. ([3]) verifying 
Dirac delta normalization. 

Comparing Eq. (fP7P with r — ► +00 to Eq. lf3Tj) with k — ► +00, one obtains: 



C, 



A4 



(38) 



The expansion of k m issued from Eq. p"?]) and Eqs. (|35I37|) . deduced from it, read in this regime: 



k — 



k n +i k n — — 
ix 



(m+|)7r q e + RV(R) p km (R) 
R 2k m R 2 Rk 2 m 

Z k {R) 



Y k (R) 
k 2 R 



A4 = l 



(39) 



where sup fc>fcn 
vanish for R — 



\Pk(R)\ has a finite limit for R — » +00 in Eq. (f3"9"P . sup k>k . \Yk(R)\ and sup fe>fe . \Z k (R)\ 
hoo. The corresponding A4 expansion is procured from Eqs. (fTTp , (|19p and (f3"9"P noticing that: 



j|(A(r)) dr 



R 



1 j t+1 (A(R))h-i(MR)) 
A(r) 



MR) 

2k 







1 - 



dr 



R Ra k (R) 
2 + fc 2 • 



(40) 



with sup fe>fe |cKfe | — > for R — ► +00. The Bessel function case, embodied in Eqs. (|24l27p , is an obvious 
particular case of Eq. (j39|) : 

^ = (m+ |)tt | a e | f3 Km (R) 
R, QifcfjiR'^ R,Hi'^ n 

it Y K (R) _ [2 Z K (R) 

where sup re>K . |/3 K (i?)| has a finite limit for R — » +00, and sup re>re . |Y^(-R)| and sup K>Km . re \Z K (R)\ vanish 
for R — ► +00 as well. 



6 Scattering wave functions in the vicinity of k = + 

We have demonstrated in Sec© that scattering wave functions converge to their open radial interval limit 
uniformly when R — > +00 for all linear momentum intervals [kd : +00 [, kd > 0. However, uniform convergence 
is generally absent on ]0 : +oo[, which is obvious when the set of positive energies wave functions always contain 
one nodeless state for R — > +00, whose energy goes to zero. It is thus necessary to devise the behavior of wave 
functions in the vicinity of k = + . For this, considered linear momenta can always be chosen smaller than 
a positive k £ momentum, which can be arbitrarily small but has to be independent of R. We will consider 
separately potentials with attractive and repulsive tails for r — > +00. 

6.1 Attractive tail 

Assuming V c < 0, it is always possible to choose R\ > Ro large enough so that the remaining Coulomb 
+ centrifugal potential is negative Vr > Ri, We can then define a fundamental ansatz solution of Eq.® 
analogously to Sec® Vr > R\. 



A„.(rl = \ J& + VQ-W±11 ,A m (r)= / ,\,„ 1, •'),//■' 



u f0 (r) = X m (r)-h ±iA ^ , u f (r) = u /o (r)[l + e(r)], (42) 

where k m is the linear momentum associated to the eigenstate of Eq.®, u(k m , r), m £ N*. 

One derives the asymptotic relation yielded by u(k m ,r) with a Neumann-Liouville power series in r~2 ; 
similarly as in Sec.©. Eq.® is used with r<j = +00 so that e(r) — > for r — > +00. Using analytical expressions 

of A m (r), A m (r) and u/ (r) of Eq.(|32")). the equality e(r) = O ( r ^) * s P rocure d f° r r — * +°° with partial 
integrations of Eq. I© . Hence : 

u(*m,r) - ^1^4™ (Jt^) 2 H Am(r) + Sm) + ° ( r ^)] ' Vr - ^' (43) 

where one has introduced the normalization A4 m , defined in Eq. (|50|) , so that Eqs. (|3QI43p are consistent, S m £ 

[0 : 2n] is the phase shift associated to u(k m ,r) and O (r~^ remains bounded for k m — > + . One has 

5 m -\ — S m — > when i? — » +00 uniformly for fc m £]0 : fe e ], m £ N*, as phase shift is a uniformly continuous 
function of k on ]0 : k e ]. 

Normalization of u(k m ,r) is defined in Eq. (fP9|) and provides, along with Eq. (j43|) . the inequality: 

-Ml m I Xmir)- 1 [sin 2 (A m (r) + S m ) + O (r^)] dr < 1 r. (44) 

App.jB]) procures an inequality from Eq. (|44|) for R sufficiently large: 



' -! - - 1 



Mr)"' dr < _ (45) 
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The k m — k m -i difference inverse is now majored. Using u(k m —i,R) — u(k rn ,R) — conditions, the fact 
that fc m _i and k m are consecutive eigenmomenta and Eq. (|43[) . one obtains: 



sin(A m _!(i?) + 8 m -i) + O(R-i) = , sin(A m (i?) + S m ) + 0{R-~-) = 
A m _x(i?) = (n - 1)tt - <y m _i + O(R-i) , A m (R) = rnr - S m + 0(iH) , n e N* 

R 

(A m (r) - A m _r(r)) dr = n + S m ^ - S m + O(R-i). (46) 

Ri 



f R dr 

u2 ,2 * ' ar 



One derives from Eq. 

(k'i - K-i) I x M "x m = ^ + - S m + 0(Br*), (47) 

where one has applied the conjugate A m (r) + A m _i(r) of A m (r) — A m _i(r) in Eq.(0B]) in order to have k% n — k 2 m _ 1 
appeared in Eq. (|4T|) . One has A m _i(r) > and 7r + <5 m _i — 6 m + 0(R~~%) > § for R large enough, so that 
Eq.(|i7|) provides: 



1 4fc„ 



k m k m — \ 7T 



R 



< — / A m (r)" 1 dr (48) 



Combining Eqs. ()45l48p . the normalization factor A/^ can be seen to be uniformly majored for k m g]0 : k e [ 
when i? — > +cxd: 

■A/fc m < M, (49) 

with M > independent of k m , k e and R. 

If V c = £ = 0, one has seen in Sec© that Eq. ([37|) is available uniformly for k s]0 : +oo], so that Eq.([49|) 
can also be utilized in this case. 



6.2 Repulsive tails 

Low-energy scattering states of a potential with a repulsive tail bear a turning point r t (k) at large distance: 



nW - *±^B+W±W (50) 

It prevents a WKB approximation similar to Eq.([42|) from being devised. However, low-energy box scattering 
states bear the property to become very close to regular Coulomb wave functions Fg n (kr) up to one constant 
when k — > + , which will be demonstrated in this section. As the set of regular Coulomb wave functions can be 
shown to be complete, using complex integration techniques and analytic properties of confluent hypcrgcometric 
functions (see Rcf. 26 for the Coulomb attractive case and App.jC]) for its extension to repulsive case), it will 
be possible to handle the contribution of low-energy box scattering states accumulating at k = + properly. 
Non-Coulombic, but repulsive tails, for which £(£+ 1) > 0, are considered in the same manner, Fe v (kr) reducing 
to Riccati-Bessel functions ji(kr) with all stated results remaining available. 

One considers < r < r t (k) fixed, r t (k) defined in Eq. (l50p and k < k t an eigenmomentum of Eq.([3]). There 
is no loss of generality therein as it is always possible to take k e small enough so that it is the case Vfe < k t . Let 
us write the equivalents of Fi v (kr) and Gt n (kr) for k — » + and r fixed, omitting factors independent of r and 
k [24]: 

F tv (kr) ~ rrMe- 7r "i^ + !(2 % /V>) , G/„(fcr) - ^^rh^k 2e+ i (2y/v^) , (V c > 0) , 

F t7] {kr) ~ (kr) i+1 , G ev (kr) ~ (kr)- e , (V c = 0), (51) 

where i 2 i+i{x) and k 2 p + i(x) are the modified Riccati-Bessel functions respectively regular at x — and x — » 
+00. 
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If r > i?o, u(r) is a linear combination of regular and irregular Coulomb wave functions: 

u(k, r) = C{k) [F tn {kr) + A(k)G ir ,{kr)} , (52) 

where C(k) 6 R is a normalization constant and A(k) £ K. a matching constant. One will show that 
A(k)Gi n (kRo) = 0(Fe v (kRo)). From Eqs. (|51l52|) . supposing the opposite would imply the Gt v (kr) component 
of u(k,r) to be dominant at r — Rq for k — » + . Then, u(k,r) oc Gg ri (kr) asymptotically for k — > + and 
i?o < r < Ri, with i?i — ► +00. This is impossible because the considered Hamiltonian in open radial interval 
bears no bound state at k — 0, so that |uo( r )| — ¥ +00 for r — > +00, with ito(r) the regular solution of Eq.([3]) 
calculated at fc = and arbitrarily normalized. 

This and the fact that -F^(fcr) and G£ V (kr) have comparable amplitudes for r > r t (k) [24j imply that norms 
and energies of box scattering states defined with Eq.Q and those defined from pure Coulomb + centrifugal 
potentials will be the same up to a value which can be made arbitrarily small Vfe G [0 : k e ] and k e small enough. 
Thus, Vk m e]0 : k e ], m <E N*, for R large enough, k e small enough and Ro < r < r t (k € ): 

\u(k m , r)\ < A u c {k' I < 2(fe OT - fc m _i), (53) 

where u(k m , r), m € N* is a box discretized scattering state of Eq.([3]), u c (k' m , r) is a box discretized scattering 
state but whose Hamiltonian is one of pure Coulomb + centrifugal parts, proportional to F£ rim (k' m r), whose 
energy differ by that of u(k m ,r) by an arbitrarily small amount, and A > is independent of k m , k' m and R 
but can be dependent of r. 

If r < Rq, u(k, r) verifies for k — > + : 

u{k,r)=V(k)(u Q (r) + o(l)), (54) 

where T>{k) is a normalization factor. As u(k,Ro) is also provided by Eq. (|S"2")) . one can match Eqs. (|52l54p at 
r = Rq, so that Eq. (|5"Tj) and the fact that A(k)Gi v (kRo) — 0(Ft v (kRo)) provides for k — > + the asymptotic 
value of V{k): 

V{k) - C(k) D rj-h-^ , (V c > 0) , 

~ C(k) D k e+1 , (V c = 0), (55) 

with Del bounded for k -> 0+. Thus, for k -> 0+: 

u(fc,r) = C(fc) ^7?-^e- 7r "KW + o(l)) , (K>0) , 
= C(fe) ^ +1 ( Uo (r)+o(l)) , (F c = 0) , 

= f(r) C(k)F in (kr), (56) 
where /(r) is bounded for k — > + from Eq. ([5T|l . Hence, the left-hand side of Eq. ([55|) also applies for r < Rq. 



7 Asymptotic treatment of bound states of potentials with attrac- 
tive Coulomb tail 

Bound states in open radial interval in our class of Hamiltonians are in finite number except if v(r) behaves as a 
Coulomb attractive potential for r — > +00, hence with V c < in Eq.([T]), which is now assumed for this section. 
As the sum of bound states in Eq. (|23[) becomes an infinite series for R — > +00, with all energies accumulating 
at zero energy, it is necessary to devise the asymptotic behavior of bound states for large principal quantum 
number to study its convergence properties. It is convenient to use n n = —ik n , with n £ N the number of 
nodes and k n the complex linear momentum of the bound eigenstate u n (r), so that n n — ► + for n — > +00. R 
is naturally demanded to be larger than Rq to have the non-local part of Eq.Q disappeared. Owing to Eq.Q, 
one defines the Coulomb turning point r[ of the bound eigenstate u n (r) for n n —* 0: 

(n) = \mvwEw±M = m + on). (57) 

ZK n K n 
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To differentiate wave functions defined in open radial interval and with box boundary conditions, one will write 
u n (r) in open radial interval and u n (r) for box wave functions. One then defines the fundamental solution 
ansatz of Eq.((3]) suitable for bound states of large principal quantum number n: 



A «W = > A„(r) = jf \ n (r>) dr' , u f (r) = A n (r)^e ±4A " «[1 + e n (r)], (58) 

with r e i > Rq, chosen large enough but independently of n, so that the square root argument in A„(r) is always 
strictly positive Vr G [r^ : r$ [ Vn > JV, N G N. Similarly to Sec. (|6.1[) . u„(r) is shown to verify the following 
asymptotic relation by way of Eq. (|58p for n — > +00: 

5»(r) = C n X n (r)-i [sin(A„(r) + 6 n ) + O (r~*)] , r < r^, (59) 

where ri^ = -3—, value introduced to avoid the divergences arising close to the turning point r[ n \ O ( r ^) 

remains bounded for n — > +00, C„ € 1 is a normalization constant and S n G [0 : 27r] plays the role of a phase 
shift for i„(r). As C~ 1: u n (r) has a finite limit for n — > +00 for r > fixed from Eq. lfST))) and as it n (0) = 
Vn G N, C~ 1 u„(r) = O(l) Vr > fixed and n — > +00. As u n (r) is normalized, one can write: 



Td fie 

■r2/„/\ j„/ 1 / -2/ / 



<(r') dr' + / u\(r') dr' < 1. (60) 



The first integral of Eq. lpjQ)) divided by has a finite limit for n — ► +00 as r<j is independent of n. Hence, the 
first integral of Eq. (f6"TJ)) is 0(C r 2 J for n — » +00. Using Eqs. (|57l58l59p and App.((B]), one obtains for the second 



dr 

Cl^ f&P f\V c \ 1 MV I 



integral of Eq. (jHO)) with n — 


+ OO: 


/ ul(r')dr> = 


Cl / Xnir)' 1 








% / Kir)- 1 

L Jr d 




0{Cl K ~ 3 ). 


Owing to boundedness of u„ 


(r) divided by C n 


r 2 


= 0(«£) , u„(r) = 



(61) 

n — ► +00 and r > fixed, Eqs. (|60|61|) provide: 

K„'j , r > fixed , n — > +00. (62) 

The following expansion is standard from quantum defect theory |32j : 

K n = 0(n _1 ) , 71 — > +00. (63) 

Eqs. (|62l63p then provide: 

u„(r) = 0(n~2) , r > fixed , n — > +00. (64) 

The case of Hamiltonian defined with box boundary conditions will be shown to be very similar for R — > +00, 
as long as n < Nb(R), with Nb(R) the number of bound states with box boundary conditions, diverging for 
R — > +00. For n < Nb(R), u n (R) — 0, so that consideration of the number of nodes of wave functions implies 
that k„+i < n n {R) < K n , with K n and k„+i defined in open radial interval, and K n (R) defined with box 
boundary conditions. Eq. (l64|) thus remains valid for box bound states, as long as n — » +00 along with R so that 
n < Nb(R) and re < R- In case ri" > -ft, Eq. lfST)) and the fact that n n (R) ~ K n for n — » +00 and n < Nb(R) 
imply that k^R = O(l) therein, re has to be replaced by R in Eq. lfBTJ)) so that using the same method as in 
open radial interval, one obtains u n (r) = O (^n n ^ for n — > +00 and n < Nb(R). Eq. ([M|) is thus still valid in 

this case. Note that n, and hence R, have to be taken large enough so that wave functions u n > (r) are very close 
for r £ [0 : Rq] Vn' > n, in order for the influence of the non-local potential of Eq.([3]) to be negligible. This 
condition is necessary, as the number of nodes of wave functions is uncorrelated with energy for the general 
non-local potential. 
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Completeness relation of bound and scattering wave functions in 
the open radial interval 



8 



The limit of Eq.j23D for R 



-oo will be shown to be, as can be expected: 



S(r,r') 



+00 f+00 - 

n=0 J ° L 



u{k,r)u(k,r') je(kr)jt(kr f ) 

TT 



dk, 



(65) 



where {u n (r)} n ^ is the set of bound states of the Hamiltonian of Eq.Q defined in open radial interval and 
{w(fc, r)} fee ]Q :+00 [ its set of scattering states, reserving notation without tilde symbol for states defined in a 
finite interval of radius R, as in Sec. 0. Scattering states are normalized with Dirac delta normalization. The 
first step is to show that series and integral in Eq. (p5|) converge V(r,r') € [0 : +oo[ 2 . Convergence of integral 
in k — » + is guaranteed by boundedness of scattering wave functions therein for fixed (r, r'). If V c < or 
t = V c = 0, this property is a direct consequence of Eq. (|4*3")) . applied with k m = k and A4 m = 1, equivalent to 
Dirac delta normalization. u(k,r) can be seen to have therein a finite limit for k — > + Vr > Ri (see Sec. (|6.1[) 
for definitions and notations). Boundedness Vr > comes forward due to u(k, 0) = equality Vfc > 0. If V c > 
or V c = 0,1 > 0, Dirac delta normalization is equivalent to have C(k) 2 (l + A(k) 2 ) = § in Eq. (|52|) . For k — ► + , 

one obtains C(k) ~ (see Sec. (|6.2p ). Boundedness of u(k, r) Vr > arises similarly as in the previous case. 
Problems of convergence, however, would occur in k = in the complex fc-plane, as the Green's function of 
Hamiltonians with Coulomb asymptotic possesses an essential singularity therein [T^ • Behavior of Eq. (p5]) at 
k — > +00 is dealt with as in Sec. ((3]). Eq. (|15p and Dirac delta normalization of u(k, r) imply that, for k — > +00: 



u{k, r)u{k, r ) Jt{kr)je{kr ) = — 

7T Z7TK 



sin[fc(r — r )] — 



V(r) + V(r') 
27rfc 



sin[fc(? 



- In] + O 



(66) 



As Eq. tfSB"]) possesses the same characteristics as Eq. (12^1) . normal convergence of the integral of Eq. (f6"5|) at 
k — > +00 with respect to (r, r') S [0 : i?t] 2 , i?t > 0, is achieved by way of a partial integration playing the role 
of Abel transformation in Sec. ||3J|. Normal convergence of the series in Eq. (|S5)) is provided by Eq. (|M|) . so that 
S(r,r') is continuous V(r, r') G [0 : +oo[ 2 . 

The difference S(r,r') — Sn(r,r') issued from Eqs. (|23l65p has to be proved to vanish for R — » +00. It is 
rewritten as: 



S(r,r')-S R (r,r') 



-E 

771=1 
+ OO 



(Nl m - l) u(k m ,r)u(k m: r')(k m - fe m _i) - B 2 m je(n m r)je(K m r')(K m - 



+ 2j[u n (r)u n (r') _ '^in(' , )'^ir l (^ , ')] + 



n=0 
+00 

E 

m= 1 



+00 " 



u(k,r)u(k,r') je{kr)jt(kr') 



dk 



u{k m ,r)u(k m ,r'){k m - k m -i) je(K m r)ji(K m r')(K m - /t m _i) 

7T 



(67) 



The series of Eq. ([23|) has been separated in two converging parts in Eq. ([67| (see Eqs. (I27I37|) and Sec. ([4])). The 
bound states u n (r) defined with box boundary conditions are finite in number, so that u n (r) = by definition 
for n > Nb(R), Nf,(R) being the number of box bound states. In case v(r) has an attractive Coulomb tail, using 
Eq. (|64p and boundedness of bound states, all functions u n (r) and u n (r), n £ N, can be majored in absolute 
value by a constant independent of R and n divided by 1 + 72.2 for r > fixed. Normal convergence of the 
bound states series with respect to R comes forward. The studied series reduces to a finite sum in other cases. 
The series thus goes to zero with R — > +00 as each of its terms obviously does. 

The two last terms of Eq. (|67p . a series and an integral, resemble a Riemann sum and its integral limit for 
infinite number of points. The fact that the sum is already a series for finite R and that the integral is improper, 
however, renders impossible the direct application of Riemann integral fundamental theorem. Demonstrating 
that the difference of the considered series and integral vanishes when R — > +00 is performed in App. (jD|) . 
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The first infinite series of Eq. ([H7|) . denoted as Si(r, r'), is eqnal to: 
+00 



(7Vf m - l) u(k m ,r)u(k m ,r')(k m - fc m _i) 



(68) 



Let k e > and m e € N* be an integer verifying fc„ le — * k e for i? — > +00. As u(k me ,r) and je(K mr r) have the 
same number of nodes Vi? > (see Sec.Q), K m£ can never be arbitrarily small for R — » +00. Eq. ([52"|) can then 
be applied to & me and n mr , so that K m5 — ► k e for i? — > +00. We will now always implicitly suppose that the 
limit R — » +00 has been effected before proceeding with fc e — > 0. 

One has Si(r, r') = S{~~(r, r') + S{ + (r, r'), whose respective general term is the same as Si(r, r'), with m < m e 
for Sf~(r,r') and to > m e for Sf (r, r'), S[ + (r,r') is demonstrated to vanish for i? — > +00 in App.((E|. It is 
thus trivially true for k e — > + . We now study Si~(r,r') by separating its Bessel and non-Bessel parts, which 
is permitted as m e is finite Vi? > 0: 



\STM\ < 



(A/^ m - l) w(/c m! r)w(fc m ,r')(fcm - fcm-l) 



m— 1 
m e _i 



^ ( B 2 Km ) jt{K m r)jt{K m r'){K m - K m -i) 



(69) 



If V c < or £ = V c = 0, one can use results of Sec. (|6.1[) . so that Eq 
i? -> +00: 



^ (A/fe m - l) u(k m ,r)u(k m ,r')(k m - fc m _i) 



m— 1 



imply for the first sum of Eq. (|59|L 

m e _i 

< J]] (A/" fe 2 m + l) |u(fc m ,r)u(fc m ,r')|(fc m - fc m _i) 

(70) 



m— 1 

m e _i 

< M ( k 



^m— 1 j 



Mfe f 



m— 1 



where M > is independent of R and fc e . Eq. ([70|) clearly vanishes for k e — > 0, result also available for the second 
sum of Eq. fMl) if t = 0. If the potential is repulsive in its asymptotic region, one is in the conditions stated in 
Sec. (|6.2[) . so that Eq. ([55|) accounts for the following inequality for k e small enough and R large enough: 



^ (Mkm _ - 1 ) u(k m ,r)u(k mi r')(k m - k m -i) 



< (Au c (k' m ,r)u c (k' m ,r') + B)(k' m - fc' TO _i), 

m— 1 

m e _ 1 

< A ^ u c (k , m ,r)u c (k' m ,r')(k' m -k' m - 1 ) + 2Bk e (71) 



m— 1 

where A > and £? > are independent of i? and fc e , and notations are the same as in Sec. (|6.2[) . Eq. ([7T1) is 
available as well for the second sum of Eq. ([69|) if t(l + 1) > 0. Let us suppose that the left-hand side of Eq. dTTj) 
does not vanish for k e — > 0, assuming also for the moment £ = 0. In Eq. ([7T|) . the sum involving u c (k' m ,r) 
functions, eigenstates of pure Coulomb + centrifugal potential, does not vanish as well. Thus, using u c (k' m ,r) 
functions instead of u n (r) and u(k, r) states in Eq. (f6"T|) , and denoting as S c (r, r') the S(r, r') function of Eq. (irj7]) 
which they generate, one obtains S c (r, r') ^ as all terms in Eq. (|67[) vanish for R — > +00 except the right-hand 
side of Eq. (|7ip . This result is contradictory with completeness of u c (k,r) and je(kr) functions, demonstrated 
independently in App. ((C|) . Indeed, with R t > fixed so that S c (r, r') is not identically equal to zero in [0 : Rt] 2 , 
one obtains Vr < Rt : 



S c (r,r') dr = 



+00 



Rt 



u c (k,r) / S c (r,r') u c (k,r) dr' dk 



je{kr) / S c (r, r') jt(kr') dr' dk = S c (r, r) - S c (r, r) = 0, 



(72) 



where all integrals inversions are allowed due to normal convergence of integral in Eq. (l65p at k — > +00 with 
respect to (r, r') S [0 : Rt} 2 . If £ > 0, assuming V c < implies that the Bessel sum of Sf~(r,r') in Eq. lGTfl) 
vanishes for k e — > 0. The remaining cases, with £ > and V c > 0, can then be treated similarly as for £ = 0. 
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Thus, Eas. (|71l72p imply that 5 z e "(r,r') goes to zero for k e -> 0, leading to Si(r,r') -» for i? -> +00. The 
equality S(r,r') — 0, V(r, r') 6 [0 : +oo[ 2 comes forward. Completeness of {u„(r)} ne N and {u(k, r)} k& ] . +oo [ is 
achieved similarly as in Eq.(|72p. with an arbitrary function f(r) bearing a Fourier-Bessel transform: 



/ f(r') S(r, r>) dr' = = V u n {r) / /(/) u n (r') dr' 



2- 



u(k,r) / /(/) u(k,r') dr' - -ji(kr) / /(/) jt(kr') dr' dk 



= y2^n{r) / /(r') Un(r') dr' + / u(fe,r) / f(r') u(k,r') dr' 

n=Q J ° J ° J ° 



(73) 



where fj(r) is the value of the Fourier-Bessel expansion of / in r. Convergence and inversions of series and inte- 
gral in Eq.(|73p are possible with f(r) decreasing sufficiently fast for r — > +00, as it implies normal convergence 

f(r+8)+f(r-8) 
2 



of series and integrals with respect to all parameters. As in the discrete case, fj(r) = lim^o 



9 Generalizations 

9.1 Coulomb divergence in r — > 

In Sec.©, the Bessel-WKB approximation could be applied for v(r) finite in r = and then seen available 
afterward for v(r) therein integrable. The important case of pure Coulomb potential, for which v(r) oc r _1 for 
r — > 0, however, has not been covered. An obvious remedy is to use Coulomb wave functions instead of Bessel 
functions for the WKB approximation in Eq.©, with v(r) is replaced by Va(r) = v(r) — in Sec.©. One can 

check that this generates an additional term in Eq.©, equal to — k k ~^°J^ \ V c . It poses no problem as it is 

0(k~ 2 ) for k — > +00, (9(r _1 ) for r — > 0, k fixed and 0(r~ 2 ) for r — > +00, fc fixed as well. The Coulomb- WKB 
approximation is obtained following the same method as in Sec.©: 

u(r) = C k F tv (kr) - ^Mlj^fcr) + O , fc - +00, 

u(r) - C k F tv {K {r)) + C k ^- , fc - +oo , r - +00, (74) 



V (r) = / v (r') dr' , A (r) = / y/k*-vo(r') dr', (75) 



<r) = C k ji(kr) + C k CL v ; ^ — #(*»0 + O ( C k -^-± ) , k -» +00, 





where notations are the same as in Eq.(fT5|). V(r) replaced by Vo(r) and sup k>kmin \a k (r)\ — > for r — > +00. As 
analytical properties of Bessel functions are used in Sec. (01, one has to express Eq. ([74]l with Bessel functions. 
Using asymptotic expansions of Fg v (kr) and derivative [24], one obtains: 

K[*(l+l) -tog(2fcr)]-\fr(r) / log 2 (fc) 

~2k 3t[kr) + °{ Ck ~ 

u(r) = ftfc (AoM + ^[^(^l)-log(2Ao(r))] + ^ + ft a£) ^ ^ f ^ ^ ^ 

where 0(k~ 2 ) is independent of r, (^(r) is alike to that utilized in Eq. ([71)l . ^(x) = j^x an d log(fc) terms 
appear in the expansion, which were absent in Eq. (|15|) . They obviously do not change convergence properties of 
series, in particular Abel transformation (see App.jAj) still provides absolute convergence to semi-convergent 
series (see Sec.©, Sec© and App.fD])). Eq. ([To]) is also no longer defined in r = due to log(2Ax) term. As 
logarithmic divergence in Eq. (|76p for r — + occurs only in wave function normalization in Eq. (|26[) . it poses 
no problem as well because it is therein multiplied by je(kr). The only other non trivial change concerns the 
expansion of k m for k m — > +00 and R — > +00 in Eq. (|39|) . which becomes: 



(m+f)7r a, + gVp(g) + RV c [V(e + 1) - log(2fc m g)] /? fcm (fl) 



In Ea.pgjl. Y k (R) and Z fc (i?) have to be multiplied by log(fc) 
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9.2 Different potential asymptotic for r — > and r — > +00 

From the conditions imposed to Hamiltonians in Sec.® and Scc. (l9.1j) . potentials bear the same centrifugal and 
Coulomb coupling constants for r — > and r — > +00: 

fe(r) - ^ ^ + — , r -»■ or r -► +00, (78) 

where /i(r) is the sum of v(r) and centrifugal potential. This is obviously a nuisance for the Coulomb part, as 
the charge of a physical system usually changes as r increases. Fixing the centrifugal potential also prevents 
the consideration of potentials such as the Poschl-Tcllcr-Ginocchio potential [27], for which £ > in general for 
r -> but £ = for r -> +00 in Eq.flTS]). 

This restriction is, however, not necessary in our method. To be consistent with Sec.©, we will continue to 
denote centrifugal and Coulomb coupling constants as (£, V c ) for r — * +00 in Eq.(|78[). but one will write them 
as (£', if r — ► 0. One can check that (£', V£) have to be used instead of (£, V c ) in Secs. (|3l4|9.ip . in order to 
treat the first turning point occurring close to r = correctly. v(r) and V(r) must be replaced accordingly by 
vo(r) and Vo(r) defined as: 

«o(r) - v(r) + W + V-fW + V _ YL , Vo (r) = f v (r') dr> . (79) 

(£, V c ) values must be employed without change in Secs. (|6l7[ ). In Sees. ([5 1 8[ ). both sets of constants (£, V c ) and 
(£', V£) have to be used. We will now see how to proceed in the two latter aforementioned sections. 

In Sec.©, (£, V c ) have to be employed until Eq. (f38]) , as one uses therein the fact that only Coulomb + 
centrifugal potentials remain after Rq. However, in order to obtain Eq. (|38j) . one compares the wave function 
of Eq.([3lj). defined with (£, V c ), with that of Eq.(fl7]l. derived with (f, VI) constants. In Eq.(|38jl. the relation 
between normalization constant Ck and A4 remains the same but (£', V£) must be utilized for the rest of Sec.©. 

In Sec.©, (£', V' c ) have to be used except for the u c (k' m ,r) and u c (k,r) functions, which were defined in 
Sec© with (£, V c ). All stated results remain valid. 



9.3 Complex potentials 

One will consider the case of non-hermitian complex potentials by way of analytic continuation of Eq. (|65p . For 
this, one introduces the following complex potentials: 

v(r) = v (r) + Xv^r) , w(r, r') = w a (r, r') + Xw^r, r') , V c = V ca + XV C1 , V' c = V' Ca + XV^, (80) 

where v ir), Ui(r), w {r,r'), wiir,r') are real potentials verifying Eqs. (|ll2[) . V co , V' Co and V Cl , V' Ci their associated 
Coulomb constants (see Eq.© and Sec. (|9.2p ) and A € C is the complex variable with which analytic continuation 
will be performed. Angular momenta and centrifugal potential are also treated with analytic continuation: 

£ = £ + X£x ,£' = £' +X£[, (81) 

where £ 6 R+, £' e R+ and £ x G E, £[ e I. < and £' must respectively verify %t{g) > -\ and K(f ) > -\ 
(see Sees. (1219.21) ). We can assume without loss of generality that the hermitian Hamiltonian defined with A = 
bears no BSEC's, as this situation will be recovered with variations of A. 

For (r, r') S [0 : +oo[ 2 fixed, Eq. (|65|) becomes a function of A, denoted as S(r, r' , A), where tilde symbols are 
suppressed for convenience and £' replaces £, in accordance with Sec. (|9.2| : 



+ OO /-+0O 

Sir,r',X) = V«„(rK(r')+ / 
n=o J o 



2, 

u(fc, r)u(fc, r') je'ikr)j£f(kr') dk 

7T 



(82) 



where A-dependence is implicitly present in all wave functions by way of Eqs. (|8QI81[) . Analytic properties 
of S*(r, r',A) depends on both analytic continuation of involved wave functions with respect to A as well as 
convergence properties of series and integral in Eq. (|82|) . 

As potentials can have complex values in Sec.©, Eq. (|6"6")) is still available, so that the integral of Eq. (|52")l 
converges for k — > +00, while convergence for k — > + of the integral poses no problem as wave functions are 
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bounded therein, due to their normalization by way of Dirac delta normalization. The series of Eg. (|82p has an 
infinite number of terms only if K(V C ) < 0. It can be handled as in Sec.®, using wave functions w n (t) defined 
as w n (t — r\V c \) — u n (r). One verifies that the use of (t,w(t)) instead of (r,u(r)) replaces V c by —1 in Eq.©. 
Thus, imaginary parts of potentials and energies are put in rest terms in Sec®, so that Eq. (164|) is still valid 
with V c £ C verifying $t(V c ) < 0. 

Existence and unicity of wave functions is provided as for the real case, using Cauchy-Lipschitz theorem and 
Fredholm theory (see Sec©), which readily defines analytic continuation in A of unnormalized wave functions. 
Wave functions singularities occur only by way of their normalization. Firstly, it is supposed that u n {r) functions 
can never become unbound or identically equal to zero with variations of A, as then they cannot be normalized 
and hence are not analytical in A therein. Secondly, analytic continuation of bound wave functions with respect 
to A implies that their normalization is effected with the square of wave functions, and not their modulus 
square. Thus, if it happens that the norm of a square-integrable bound state is equal to zero, the state cannot 
be normalized, i.e. it is an exceptional point [28 . They appear with complex potentials only and are obtained 
by making two different bound or resonant states of a real potential becoming degenerate by way of increase 
of the potential imaginary part [28]. This obviously occurs only for a discrete set of complex A's. Another 
source of wave functions poles is induced by spectral singularities [29] ■ Contrary to real potentials, complex 
potentials can bear resonant states of real energy E > 0, for which S~ = in Eq. ([5D|) . This is incompatible 
with the requirement that S + S~ = \ (see Sec.©), required by Dirac delta normalization, so that they generate 
a singularity at k 2 = E. Contrary to exceptional points, they can be handled by deformation of the A:-real axis 
in the complex plane, where application of Cauchy theorem have them appeared as poles, similar to resonant 
states: 



S(r,r',X) = V u n (r)u n (r') + V u$(r)u$(r') + 



71=0 



£ 

m— 1 



L- 



2 ~ 

u(k,r)u(k,r') ji'(kr)jt(kr') dk 

IT 



(83) 



where M 6 N* 



r) is a properly normalized spectral singularity lying under the complex contour L + , which 



can be the real axis except in the vicinity of spectral singularities, where it consists in small half-circles in the 
lower complex plane encompassing Um\r) states. This procedure allows to extend the domain of analyticity 
of S(r, r' , A), and in fact generates full Berggren completeness relation [14] if one allows «m(r) to be resonant 
as well (see also Ref.[30] for a study of short-range complex potentials spectrum completeness properties with 
Berggren completeness relation). If (r) states move to the upper part of the complex plane for A close to 
that generating spectral singularities, they become bound states so that L + can be replaced by the real axis, 
recovering Eg. ([82]) . but with u$(r) states being added to the series of Eq. ([82]) . For the A value generating 
spectral singularities, the only possibility to generate a complete set of states is to leave L + in the complex 
plane and to use Berggren completeness relation. BSEC's can be treated in a similar manner (see Ref.[3T] for 
a study of BSEC's induced by non-local potentials.). 

Hence, series and integrals Eg. ([52"]) converge normally with respect to (r, r') G [0 : Rt] 2 , Rt > and A for all 
A's preventing singularities from appearing in series and integral. Handling of singularities crossing the real axis 
is effected with aforementioned contour deformation. Thus, the analyticity domains of S(r, r' , A) and involved 
complex wave functions are the same. If A G K, one is in the hcrmitian case, so that S(r, r' , A) can be treated 
with Eg. ([55]) . Hence, it is identically equal to zero in a small interval [Aj,A/] containing A = 0, chosen so 
that no bound state of energy zero occurs in this interval. S*(r, r' , A) = immediately follows in its domain of 
analyticity. Completeness of wave functions generated by a complex potential is then procured as in Sec.©, as 
the set of Riccati-Bessel functions for £' 6 C, 3F£(£') > is complete (see App.JCj). 



10 Conclusion 

Restricting Hamiltonians to vanish identically after a finite radius Rq, except for a pure Coulomb potential part, 
has allowed to demonstrate the completeness of their eigenstates by way of analytical properties of Bessel and 
Coulomb wave functions. It has also allowed non-local potentials to be included, as their localization before i?o 
renders their effect negligible for k — » +oo. The subtraction of two completeness relations of wave functions 
defined with box boundary conditions, one involving the studied eigenstates and the other generated by Bessel 
functions, has made possible to approach the problem with elementary methods, without recurring to Lebesgue 
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measure theory. The fact that the obtained series converges has a strong physical significance: completeness of 
eigenstates is directly related to their asymptotic convergence to sine and cosine functions for k — > +00, whose 
completeness properties are those of Fourier series and are obvious from Dirichlet kernel properties. The singular 
point at k — > 0, due to the presence of the Coulomb potential at r — > +00, possesses the same characteristics as 
that of pure Coulomb Hamiltonian, whose wave functions are analytical expressions of confluent hypergeometric 
functions. Their completeness relation can thus be studied independently with complex function theory. The 
open radial interval case can then be procured with the box radius sent to infinity, where all limiting processes 
can be handled properly due to asymptotic properties of eigenstates for k — > and k — > +00. 

Far from the drastic approximations usually effected on Hamiltonians to derive simple demonstrations of 
completeness, the set of considered Hamiltonians covers almost all physical situations assuming spherical sym- 
metry and pure Coulomb asymptotic, the difference between a potential decreasing quickly to zero or vanishing 
after a finite radius being nothing but academic for practical purposes. It is hoped that this demonstration, 
derived with full rigor but accessible to non-specialists in operator theory, will provide better insight to the 
difficult problem of Hamiltonians completeness relations of eigenstates. 
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A Abel transformation of semi- convergent Fourier series 

The partial sums defined by the terms of Eqs. (128[) generating non-absolutely convergent series can be noticed 
to be of the form: 

y /(x) sin(mx+a) , (ai) 



where 



or 



R ' \ 2R / 2rmrR ' V ' 



R ' V 2R J 2mirR y ' 

In these definitions, M e N*, x S [0 : 2ir] in Eq. (|A2p while x 6 [-n : ir] in Eq. (|A3|) . f(x) is differentiable due 
to Eq. (TT6]) and f(x) = if sin(ma; + a) = 0. It is then convenient to use the following function g{x) and partial 
sum Sm(%)' 

= f(xy {x+a) {x ^ {0j 27r}) (s) = ie i af{x) {x e {0j 2?r})j (A4) 

1 — e lx 

^ pimx 

S M {x) = g(x)e- ix (l-e ix )Y, , (AS) 

rn—l 

where g{x) is continuous on its domain of definition and Eq. (|Al[) is equal to the imaginary part of Sm(x) 
VM G N*. One can apply Abel transformation to Eq. l|A"5|) for M e W: 



M /I 1 \ 1 — p iMx 

S M ( X ) = g(x) E (1 - e-) - - + fAx)^^ (A6) 
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Letting M — > +00 in Eq. (|A6[) , one obtains the value of the its limit series S(x): 

+ OO , s 

^)^WEd-0 --^ (A7) 

m— 1 v ' 

S(x) is normally convergent in Eq. (|A7[) with respect to x and the series defined in Eq. (|Aip with M — > +00 is 
equal to its imaginary part. As a consequence, Eq. (|28p can be rewritten so as to generate a normally convergent 
series V(r,r') G [0 : R} 2 . 



B Equivalent of norms of wave functions for V c < and k —> 



In Sec. ()6.ip and Sec. (O, one has to devise the equivalent of the integral X for 77 — > +00, representing a part of 
the norm of a low-energy wave function u(r) in an attractive Coulomb field: 



1 - 



u(r) oc A(r) * sin(A(r) + 5) + O fr *J 
A _1 (r) [sin 2 (A(r) + 5) + O (r~^\ 



dr. 



(Bl) 



Eq. (|Bl[) becomes after partial integration 
2 



' 1 ' \- l {r) dr - ~[sin(2A(r) + 26) \- 2 (r)] r r ( 



+ 



2 

1 /■''' 



d 



m(2A(r) + 25)—(\(rT 2 ) r , =r dr + 



O 



X-\r) 



dr. 



(B2) 



One will show that Eq. (|B2|) is equivalent to its first integral if rt — > +00. 

Let us first consider u(r) to be a scattering state of linear momentum k defined with box boundary conditions 
(see Sec QG-ll) 1 ). In this case, one has: 



/■, = R , r f = R, A(r) = y k 2 + ^ - , A(r) = J* A(r') dr' 

where i?o is fixed and i? — >• +00. Eq. (|B3[) implies for 77 — ► +cx) and fc E]0 : fcj, /c e > 0: 



A- 2 (r 7 ) - 0(ry) , ^ (A- 2 (r)) r=r/ - 0(1) , J \-\r) dr = O 
If u{r) is a bound state of linear momentum k n (see SecQ), one has: 



'■f- 



ri = r d ,r f 



- „(r0 



A(r) = W^ 



4 , A(r) = / A(r') dr' 



(B3) 



(B4) 



(B5) 



where is fixed, n n — —ik n and = + 0(1) is half the value of the turning point r[ of u(r), defined 
in Eq. (|57|) for K n —> 0. Owing to Eq. (|B5|) . one obtains, with K n — ► 0, and thus r/ — > +00: 



A(r/) ~ , A- 2 (r ; ) ~ K " 2 , - (A- 2 (r)) r=r/ = 0(1) 



A (r) dr = K r 



dr 

t 



i{t + 1) 
t 2 



<ft = 0(«" 3 ). 



(B6) 



For situations embedded in Eqs. (|B4IB6[) . the following properties hold for r/ — * +00: 



\-\r) dr^+^, \- 2 (r f ) = o ( / A^r) dr ) , - (A- 2 (r)) r=r/ = ( / A _1 (r) dr 



rr t \-\r) 



dr = o 



A _1 (r) dr 



(B7) 
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Standard theorems can then be used with Eqs. (|B2IB7[) to provide the result: 



A _1 (r) sin 2 (A(r) +6) + (r~*) dr ~ i J ' A" 1 (r) dr , r f -> +00. (B8) 

C Completeness relation of Coulomb wave functions in a repulsive 
field 

In our demonstration, the completeness relation of Coulomb wave functions generated by a pure repulsive 
Coulomb + centrifugal potential is required. The attractive Coulomb + centrifugal potential with I S N has 
been treated in Ref. [26] . but the repulsive case was not mentioned. However, due to the very similar structure 
of attractive and repulsive Coulomb wave functions, the demonstration of Ref. [26] just has to be mildly changed 
in order for the repulsive potential to be taken into account, as well as for feC, $l(£) > — A. 
Following the method used in Ref. [26] . one defines the integral: 



2 f K 

J(K;r,r') = - F lv (kr)Fe v (kr') dk (CI) 
t Jo 

where K > 0, r > 0, r' > 0, 77 is the Sommerfeld parameter equal to V c > and Coulomb wave functions 
are properly normalized to the Dirac delta normalization. One demands also that r' > r, which can be clearly 
effected without loss of generality. For K — > +00, this integral will converge weakly to S(r — r'), as no Coulomb 
bound state appears in the repulsive case. Using the expression of Fe n (kr) in terms of confluent hypergeometric 
functions [24], denoted as $(a, b, z) and ^f(a,b,z) for regular and irregular functions respectively, Eq. ljClj) 
becomes [26] : 

J(K; r, r') = (±>j ^£+7)* f ^ + 1 + l7 ^ £ + 1 ~ iv)k 2i+2 e ik ^e~^ 

x <P(£+l + ir],2£+2,-2ikr)<P(£+l + ir],2£+2,-2ikr') dk (C2) 
— ^ ( 4rr 'y +1 f K k 2i+2 etkr ^ e +1 + l11j2 £ + 2, -2ikr) 



2-k J T{2£ + 2) J 
x f e -Ml+i)+ikr'T(£ + 1 + i v )^(£ + l + ir ]l 2£ + 2, -2ikr') 

+ e l ^ +1) - lkr 'T(£+l-ir,)m{£+l-iri,2£ + 2,-2ikr'e™)} dk (C3) 

- (s) *w L ^""^ (' + 1 + 

x T(£+l + ir))^(£+l + ir],2£+2,-2ikr') dk. (C4) 

Eq. ([C3p is obtained making use of the definition of ^(a, b, z) as a linear combination of $(a, b, z) functions for 
2£ $ N, the case of 2£ e N being accounted for by a standard limiting process [24]. Negative fc's in Eq. (]C4p 
are treated as \k\e l7V [55], as one will consider in the following complex fc's belonging to the upper part of the 
complex plane only. Note that k 2l+2 is no longer an even function of k if £ N. 

The integral of Eq. (]C4[) runs from —K to K, so that one can use complex integration with a half-circle 
contour of radius K in the upper part of the complex plane to which the real axis is added. The irregular point 
at k = will be treated using a small half-circle of radius k e in the upper part of the complex plane which will 
be shown to vanish for k e — > 0, similarly to Ref. [13], e ife ' r+r ■* $(l + 1 + ir), 2£ + 2, —2ikr) has a finite limit for 
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k — > 0, directly determined from $>(£ + 1 + irj, 2£ + 2, —2ikr) power series expansion 



3l fc(r+r') + ! + 2 £ + 2, -2ifcr) = 



' + ''> Shroff- 2 *'" 

n— 

(rV c - 2ifcr(f +!))••■ (rF c - 2iAx(£ + n)) 



a ik(r+r') 



n=l 



n! {21 + 2) n 



E 



n! (2^ + 2) 



(C5) 



where (a)„ = a • • • (a+n— 1) is the Pochhammer symbol and the series limit at k = is effected summing the limit 
of all series terms at k = 0, which is permitted due to normal convergence of the series in the unit circle of the k- 
complex plane. The remaining term of the integrand of Eq. (IC4l) is treated expressing fy(£ + 1 + irj, 2£ + 2, —2ikr') 
as a linear combination of $'s if 21 $ N or with its logarithmic expansion if 21 g N [2H [55]. The integrand of 
Eq. (|C4p is thus bounded in a complex neighborhood of k = 0, so that the contribution of the mentioned small 
half-circle integral vanishes for k e — > 0. 

Poles of Eq. (|C4p are generated only by T(i + 1 + irj), similarly to Ref. 26 . They are immediately seen to 
bear negative imaginary part, and thus no pole arises inside the used complex contour, contrary to the situation 
of Ref. [2"6"] where they correspond to Coulomb bound states. 

For K — ► +oo, the integral on the half-circle of the complex contour can be calculated with expansions of 
confluent hypergeometric functions for fc| — > +oo [2"4l I26j : 



^(£+l + irj,2£ + 2, -2ikr') 



{^ikr'y 1 - 1 -^ 1 



[£-ir])(£ + l + ir]) 



{-2ikr 



.VMk) 



-2ikr' 



Ak- 



T(£ + 1 + in) 
T(2£ + 2) 



$(£+l+ir},2£ + 2,-2ikr) = e 



- ^^X^ikr)-'- 1 



+ e- 2ikr (-2ikrY 



2k 

1 — % 

2k 

■ Klog(fc) 
' l ^2k— + 



L + o 

+ Bk~ 1 ■+ 
Ck^ + O 



+ 0(k- 2 ) 

log(fc) 2 
k 2 

O 



(C6) 



log(fc) : 



log(fc) 



(C7) 



where "±" is the sign of — 3?(fc) and A, B and C are complex constants independent of k. 

Let us denote as C(K;r,r') the aforementioned integral, so that k therein runs on the complex contour of 
the half-circle, i.e. k = lie 10 with 9 £ [0 : it}. The integrand of C(K;r,r') is naturally the same as in Eq. (|C4[) . 
Inserting Eqs. (|C6IC7[) in Eq. (|C4[) . the expression of C(K ; r, r') for K — > +oo comes forward: 



C(A>,r') = 



s -tir(^+l) 

2^ 
2n 



2 i^e 4Keie(r ' +r) 



,y c log(^) + ^_ 1 



A" 



gW _ .I/ c log(AV*) 



- -rEK- 1 
log(A-) 2 



d0 



— C iKe iKei "^-^ [e^ + FK- 1 ] d9 + O 
2 71 " Jo 

e -x(r+r') ^(7,.^ + !)) + s i n (A-(r + r') - n(£ + 1)) 
7r(r + r') 

sin(A(/-r)) F 

— 5„, +o(l), 

7r(r — r) 2 



(C8) 



(C9) 



where _E and F are constants independent of K and and 8 rr i is a Kronecker delta. Care was taken in handling 
powers of complex numbers when deriving Eq. (|C8j) . using the fact that all k's have positive imaginary parts. 
The r' > r condition is necessary for the exponential depending on r' — r in Eq. (|C8|) to remain smaller than 
one in modulus. Eq. (|C9p implies that C(K; r, r') — ► —S(r — r') weakly when K — » +oo by way of its second and 
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third term, its first term being proportional to S(r + r'), vanishing as r + r' > 0, while all the others are shown 
to go to zero with integration by parts. Cauchy theorem then provides the result: 

9 r+oo 

- F er ,(kr)F eri (kr') dk = 5(r-r'). (CIO) 
* Jo 

Due to the appearance of Dirichlct kernel in the demonstration, Fourier-Coulomb transform obeys the same 
convergence behavior as Fourier transform, i.e. the Coulomb-Fourier expansion of a function f(r) is equal to 

can here be equal to zero, the completeness of spherical Riccati-Bessel functions in 
open radial interval has been demonstrated as well. 



D Convergence of series to improper integral via Riemann integral 
fundamental theorem 

The expression of interest in Eq. (pT)) is the following: 



+00 
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(Dl) 



It is conveniently rewritten as: 
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771 — 1 

^ ^ [J{k"m){km ~ km— l) ~~ J (^m)i,^m ^m — l)] ) 



(D2) 



where U(k) and J(k) are respectively u(k, r)u(k, r') and ^ji{kr)jt(kr'). Infinite series in Eq. (|D2|) are convergent 
due to Eqs. (|22I24I66|) . Let us consider K > and R > and define M K E N* so that k Ml< -» K for R -» +00. 
Reasoning with km k as with n mt in Sec®, one shows that km k — * K for R — > +00. 
We will focus first on the rests of the series of Eq. (|D2|) with respect to Mk ■ 



+ OC 



Ru = ^ [U(k m ) ~ J(k m )](k m - k m -i), 

171 — Mk 

+ OO 

Rj = ^ [J(km){km ~ &m-l) — J{^m){^"m ~ ^"m— l)] 



(D3) 
(D4) 
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It is convenient to transform the series of Eqs. (|D3ID4[) by way of an Abel transformation (see App.((X|) so that 
they bear absolute convergence. We first study Eq. (|D3[) . Using Eq. ([66)l . the semi-convergent part of Eq. (|D3[) 
can be written as the imaginary part of the sum of two terms of the following form: 



+ 00 



m=M K ^ m " m+1 ' n=M K 

where x and f{x) are equal respectively to r±r' and T[V(r) ± V(r')]/27r, with the last term multiplied by e~ llT! 



Fu{x)= A ^ X nT--ir—) > A ™(X)= E f( X ) eiknX ( k n ~ K-l), 



(D5) 
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if x = 



r'. A m (x) in Eq. (|D5[) has an simple asymptotic expansion for R — > +00, deduced from Eq. 



f{xy k - x 



=M K 
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2i?ri7r 
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DO,' 

where sup fc>i ^ | ctfe | — ► for R — » +00 and the fact that Y^ n {R) = 0((k„ 
been used. The part of the sum in Eq. (|D6|) depending on can be handled with Abel transformation as well. 
From Eq. (|D6p . A m (x) is bounded for k m — > +00 and 1 £ R, as is finite for x — > 0. As a consequence, Abel 
transformation performed on Eq. (|D3p provides the inequalities: 



log(-R)) for R — > +00 has 



< A> E 
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+ 00 
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(D7) 



where the property that k m > fc m _i Vm G N* has been used and Aq, A are positive constants independent of 
K and R. Using Eqs.j22]) and (J39]) , one can rewrite Eq. (jD4|) as: 
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2irR 

-\-oo 
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ttl—Mk 
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fc"m— l) 



(D8) 



where sup K>if |a K (i?)| — > for i? — > +00. Abel transformation of the semi-convergent term of i?j in Eq. (|D8[) 
can be accounted for in the same manner as for Ru- Consequently, \Ru \ + \Rj\ can be majored by 4>, where 
A is a positive constant independent of K and R. We have demonstrated that Eq. (|D2[) can be majored by the 
positive constant P: 



M K 



I [U(k) - J(k)] [U(k m )(k m - k m -i) - J{n m )(n m - K m _i)] 

J ° m=l 
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(D9) 



Let us take e > 0. The two last terms of Eq. (|D9[) can be made smaller than e for K large enough, independently 
of R. The remaining term of Eq. (|D9p is the difference between a Riemann sum and its integral limit on the 
finite interval [0 : K] for R — ► +00. Thus, it can made smaller than e for R large enough. Eq. (|Dip has thus 
been shown to go to zero when R — ► +00. 



E Limit of series S^ + (r, r') for R — > +00 

One considers the series Sf + (r,r') defined in Sec.©: 



Si + (r, r') = [ ( M L ~ !) u(k m ,r)u(k m: r')(k m - k n 



B k h(K m r)je(K m r')(K m - 



(El) 
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where u(k m ,r) is defined by way of Eq. (f6"5)) in Sec.©, ji(K m r) was introduced with Eq. (j2T3)) in Secfg]), m e 6 N* 
is defined so that k mt — > k e and K m<! — > k e for R — > +00, fe e being a positive linear momentum (see Sec.®) and 
B Km , Mk m are respectively defined in Eq. (|2"Tj) and Ea.([3"D)). 

Due to boundedness of scattering wave functions and absolute convergence of series (see Eq. (|2"T)) ). \Sf (r, r')\ 
verifies the inequality: 



+00 +00 „ 

B L - 

IT 



\Sf + (r,r')\ < M u WL ~ !K fc ™ ~ k m-i) + ^ 



(Km ~ «m-l), (E2) 



where M u and Mj are two positive constants independent of the box radius R. From Eqs. (|39l41[) . A4 m and 
B n m verify: 



a 1 1 ^ M M B ) 



L- 2 



9 2 



where fc m > fc e , K m > fe e , and M^j-(R), Mb(R) are two positive constants independent of k m and K m and 
vanishing for R — > +00. Eq. (|E2p implies for £ ; £+ (r, r') and i? sufficiently large: 

+00 , , +00 
\St + (r,r>)\ < M U M N (R) £ * m ~ * m ~ l + M B (R) ]T ^^zi 

m— m e m—m e " L 
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< M U M M (R) h , m ' 1 +M j M B (R) ]T 

+ OC , +00 , s 

< M u Mf/(R) £ +M i M fl (i?) E (_-_) 

rr)=rn _ x ' 7Tl—Tn e 



M u M M (R) + M 3 M B (R) 



< 2 " Ay v ' J flV " ; , (E4) 
where one has used the properties k m > fe m _i and K m > K m -\. Hence, Sj (r, r') — » for i? — > +00 V fc e > 0. 
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